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Abstract. - The dynamic critical exponent z is determined from numerical simulations for 
the three-dimensional XY model subject to two types of dynamics, i.e., relaxational dynamics 
and resistively shunted junction (RSJ) dynamics, as well as for two different treatments of the 
boundary, i.e., periodic boundary condition (PBC) and fluctuating twist boundary condition 
(FTBC). In case of relaxational dynamics, finite size scaling at the critical temperature gives 
z ~ 2 for PBC and 1.5 for FTBC, while for RSJ dynamics z ~ 1.5 is obtained in both cases. 
The results are discussed in the context of superfluid/superconductors and vortex dynamics, 
and are compared with what have been found for other related models. 



A neutral superfluid like 4 He and a superconductor in the limit of large London penetration 
depth can be characterized by a complex order parameter, and the XY model can be viewed as 
a discretized version of this type of systems Q . All these systems are expected to belong to the 
same universality class for the thermodynamic critical properties of the phase transition. An 
interesting feature of these models is the presence of thermally generated topological defects. 
In two dimensions (2D) the topological defects take the form of vortices and give rise to the 
Kostcrlitz-Thoulcss transition 

§ 1- Also in 3D thermally generated vortex loops are present 
at the transition and it has been argued that the critical properties, both the static and the 
dynamic, can be associated with these vortex loops || |J, ||. The low-temperature phase in 
the 3D case consists of closed vortex loops of finite extent whereas in the high-temperature 
phase the loops can disintegrate [|[ |J [|. 

In the present Letter we investigate the dynamic critical properties for two simple cases 
when the static properties are given by the 3D XY model. One connection between the 
vortex loops and the dynamical properties is through the 2tt phase slip across the system 
occurring when a vortex loop expands so much that it leaves the system. The connection is 
most easily phrased in case of a superconductor: the rate at which the vortex loops expands 
and leaves the system, when driven by a dc current, is proportional to the voltage across the 
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sample. Consequently the vortex loops are connected to the resistance which for a system in 
equilibrium can be obtained by the fluctuation-dissipation theorem Q. 

Similarly to the static case one expects universality also for the dynamic critical properties 
in the sense that the critical dynamics does not depend on the details but rather on general 
characteristics like conservation laws, spatial dimensions, and the static critical properties 
The motion of the vortex loops are associated with a conservation law since the vorticity of a 
fixed area can only change by vortex segments leaving or entering the area. This conservation 
law restricts the motion of the vortex loops and consequently one may expect that the longest 
relaxation time of the system can be associated with the vortex loops. In the dynamic 
universality classes defined by Hohenberg and Halperin (tJ the dynamics of a 3D superfluid 
belongs to model F characterized by the dynamic critical exponent z = 3/2[|). A model with 
purely relaxational dynamics on the other hand belongs to model A with z«2 |. 

The two cases we study are relaxational dynamics which does not have local current 
conservation, and the resistively shunted junction (RSJ) dynamics which has local current 
conservation. The Hamiltonian for the 3D XY model on an L x L x L cubic lattice can be 
expressed as 



where the sum is over all nearest-neighbor pairs, — 6j — • A is the difference in the spin 
direction between the neighboring sites i and j, and J is the coupling strength. The twist 
variable A = (A x , A y , A z ) is a vector such that LA X measures the average rotation of the 
spin direction when going from one boundary surface to the opposite in the x direction and 
similarly for the other directions, is the unit vector from site i to the nearest-neighbor site 
j (the lattice spacing is taken to be unity), and 9i is a phase angle associated with each site 
i measured with respect to the local spin-direction associated with a uniform twist A across 
the sample. We use the boundary conditions 6{ = Oi+Lz = Qi+Ly = ®i+Li- 

In the superfluid/superconductor analogy of the XY model, 6*; — 9j — r,y • A is the total 
(gauge invariant) phase difference of the order parameter and the twist A may be thought of 
as the contribution to the gauge-invariant phase from a spatially uniform vector potential. 

The usual periodic boundary conditions (PBC) for the XY model correspond to periodic 
spin directions and to A = 0. However this imposes an unphysical restriction on the topological 
defects (the vortex loops): the original defect-free state is not regained with PBC, when a defect 
is created in a defect-free state and then annihilated across the boundary [^) . The additional 
degrees of freedom introduced by A, on the other hand, ensure that the energy associated 
with a given configuration of topological defects is unique ||. The more physical boundary 
condition which includes the fluctuations of the twist is termed the fluctuating twist boundary 
condition (FTBC) 

In the RSJ case the total current from i to j is the sum of the supercurrent, the normal 
resistive current, and a thermal noise current: 



where i c = 2eJ/h is the critical current of a single junction, Vij is the potential difference across 
the junction, r is the shunt resistance. The current conservation law at each site, together 
with the Josephson relation d(8i — 0j — Vij ■ A)/dt = 2eVij/h, allows us to write the equations 
of motion in the form 




(1) 



i c sm(6i - 6j - • A) + — + rjij, 



(2) 




(3) 
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where the primed summation is over six nearest neighbors of j, Gij is the lattice Green function 
on the cubic lattice. For convenience we from now on use units such that i c = J = r = h/2e = 
1. The remaining dynamical equation for A is obtained from the local current conservation 
together with the global current conservation condition that no currents pass through the 
boundaries pj: 

A = - rA SA + " A (4) 

with Ta = 1/L 3 . In order to ensure the correct thermal equilibrium the noise correlations obey 
the relations: (r]ij(t)) — 0, (0)) = 2T(SikSji — Su8jk)S(t), and correspondingly for the 

three components of rj^: (r]/^ m (t)) = 0, (r)A m (t)VA n (0)} = (2T/ L 3 )S mn 5(t) with m,n = x, y, z. 

The RSJ equations defined in this way incorporates local current conservation and the 
boundary conditions are chosen such that the vorticity for each of the six sides of the cubic 
lattice is zero at any instant and that there is no current flow through them. The RSJ equations 
are usually phrased in the superconducting language, however, they also apply to a neutral 
superfluid; the RSJ equations in this case correspond to a constant mass density and local 
conservation of mass current. 

We also note that eq. (^), which is related to the resistance and hence to the vortex loops, 
has a relaxational form. However, one should note that the relaxational constant Ta = L~ 3 
is unusual since it vanishes with the size of the system and that eq. (|]) by itself can be viewed 
as a global current conservation law expressing that the average total current of the system 
vanishes at each instant. 

How important is strict local current conservation for the critical dynamics? We investigate 
this by comparing the results from the RSJ dynamics to relaxational dynamics, where eq. (|^) 
is replaced by the purely relaxational form: 

Mi(t) OH 

-^r = - T w i +r]iit) (5) 

with (r)i(t)) = and (r]i(t)rij (0)) = 2TSijS(t) (we have set T = 1). Thus the dynamics is 
in this case given by the two relaxational equations (Q) and (|^). Superficially one might 
have guessed that this relaxational dynamics should belong to model A|j. However, as will 
be shown below, the z value obtained from size scaling at T c is not compatible with this 
expectation. This suggests that the global conservation law reflected in the size dependent 
relaxation constant in eq. (^) is enough to slow down the critical dynamics. 

The dynamical equations are integrated numerically using the second order algorithm in 
ref. Jll| with a discrete time step At = 0.05 for RSJ and At — 0.05 and 0.01 for relaxational 
dynamics, using lattice sizes up to L — 32 and 24, respectively. 

The resistance R is related to the equilibrium fluctuations of A(i) by the fluctuation- 
dissipation theorem JhJ: 

^=~<[A m (e)-A m (0)] 2 ), 

where G is a large enough time interval (9 = 2000 in the present simulation). Near the second 
order phase transition the intensive quantity LR(T, L) obeys the scaling relation: 

LR(T, L) = L-^-VplL^^T-Tc)]. (6) 

Thus if we plot the ratio ]n[R(T,L)/R(T,L')]/ln(L/L') as a function of T for different pairs 
(L,L r ) then the curves should cross at (T c , — z) |13]. The inset in fig. 1 shows the RSJ result 
for the pairs (L, L 1 ) =(4, 8), (8, 16), and (4, 16). The crossing point gives T c = 2.20, which is 
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Fig. 1. - The scaling curve for resistance 7? in the RSJ case. The parameters used are z — 1.46 and 
T c = 2.20, determined by the intersection construction described in the text and shown in the inset, 
together with the value v — 0.67 expected for the 3D XY model. From these scaling constructions 
we obtain the estimate z = 1.46 ± 0.06. 



very close to the true critical temperature T c w 2.202 for the 3D XY model fl4|j , and z = 1.46. 
Figure 1 confirms this determination by aid of the full scaling relation eq. (H) using T c and 
z obtained above together with v = 0.67 for the 3D XY model {v « 0.671 ± 0.001 |||). 
As seen a very good scaling collapse is obtained. We have estimated the precision in the 
determination by treating z and T c as free parameters in the full scaling relation with the 
result z — 1.46 ± 0.06. One reason for treating T c as a free variable is that in principle the 
finite time step At in the integration introduces an uncertainty in T c [ [To|| . Figure 2 gives the 
corresponding result for relaxational dynamics. In this case the integration turns out to be 
more sensitive to the choice of At. In order to handle this we calculate R at T c — 2.20 where 
R oc L~ z [see eq. with At = 0.05 and 0.01, extrapolate linearly to At = 0, and obtain 
zwl.5. 

Our conclusion is that RSJ dynamics and relaxational dynamics have the same size scaling 
of the resistance at the critical temperature and that the z value obtained from this scaling 
is consistent with z m 1.5. This implies that the constraint imposed by the local current 
conservation is less crucial than one might have thought. We suggest the following hand- waving 
explanation: Equation (^) describes individual spins relaxing towards a state with a given 
value of A. Each fixed configuration of vortex loops correspond to a twist A. Since A has 
a slow relaxation governed by eq. (Q) this suggests that the dynamics is compatible with a 
situation where the change of a vortex loop configuration is slow compared to the relaxation 
of the individual spins. From this perspective RSJ dynamics and the relaxational dynamics 
described above are just two alternative ways of imposing a slow dynamics on the vortex loops. 

Next we consider the case when A = which corresponds to the standard PBC imposed 
on the spins. Relaxational dynamics is in this case given by eq. Q with A = 0. This is 
compatible with the situation when the spins are relaxing directly towards the global ground 
state with A = 0. In this case we cannot use eq. (ph to find z (because A = 0). Instead we 
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Fig. 2. - Determination of the dynamical critical exponent z for the relaxational dynamics with FTBC 
from resistance scaling R oc L~ z at T c = 2.20. The data obtained for L = 4, 8, 10, 12, and 16 and 
the two integration time steps At = 0.05 and 0.01 were linearly extrapolated to At = 0. The value 
z = 1.47 was obtained from a least square fit. 



resort to the following size scaling relation valid at T e [fl6| 



where G(t) — (F '(t) F '(0)) / L 3 is the time-correlation function, with F(t) = sin(0, — 9j) 

and the sum is over all links in one direction EQ] . In the superconductor analogy this is related 
to the supercurrent correlations [lQ|, Figure 3 shows that a good scaling is obtained for z = 2. 
From this we conclude that relaxational dynamics with standard PBC has z ~ 2 consistent 
with the model A universality class jj], ^2|. Our suggested interpretation is that, when the 
average twist A is not one of the dynamical variables, the global current constraint reflected 
in eq. (^) is no longer present and the relaxational dynamics becomes of standard relaxational 
type. We compare this to the RSJ dynamics for the same case, i.e., standard PBC for which 
A = 0. Figure 4 shows that a good scaling collapse is obtained for z — 1.5, which is the same 
as for FTBC. This suggests that local current conservation is a sufficient but not a necessary 
condition for imposing the slow vortex loop dynamics. 

The suggestion, that the dynamics of the vortex loops can be associated with an exponent 
z rs 1.5, can be further substantiated in the following way: The Hamiltonian for the XY 
model with FTBC is dual to the lattice vortex loop model with PBC|J. This means that 
relaxational dynamics for the XY model with FTBC corresponds to relaxational dynamics 
for the lattice vortex loop model with PBC. In the latter model the only degrees of freedom 
are the vortex loops and it has been shown that z w 1.5 is obtained from size scaling of the 
resistance for this model|^, |l3[ ]. From this perspective it is tempting to conclude that z w 1.5 
can be associated with pure relaxational dynamics for the vortex loops. In addition one may 
note that an exponent z w 1.5 is also consistent with z = 1.44 in ref. || obtained from a 
theoretical treatment of vortex loops. 

Finally we note that in ref. Q the value z = 1.5 ±0.5 was determined for the RSJ model in 



G(t) = -h(tL- z ), 



(7) 
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Fig. 3. - Determination of the z for relaxational dynamics with PBC from the scaling function eq. (|7|) 
LG(t) = h(tL~ z ): A good scaling collapse is obtained for z = 2. The inset shows that no scaling 
collapse is obtained for z = 1.5. 




Fig. 4. - Determination of the dynamical critical exponent z for the RSJ case with PBC from the 
scaling function eq. ([?]) LG(t) = h(tL~ z ): A good scaling collapse is obtained for z = 1.5. The inset 
shows that no scaling collapse is obtained for z — 2. 



the presence of external currents and that in ref. jl9| z = 1.38±0.05 was found from simulations 
of a version of the XY model with spin dynamics which is an alternative dynamics consistent 
with superfluids. For the lattice vortex loop model with Monte Carlo dynamics z — 1.45 ±0.05 
was obtained in ref. |l3| and z = 1.51 ± 0.03 in ref. [pj using the same method as in the inset 
of fig. 1. 
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This leaves us with the following two main alternatives: z determined from size scaling 
for 3D XY model with relaxational dynamics and FTBC, the lattice vortex loop model with 
relaxational dynamics and PBC, as well as the 3D XY model with RSJ dynamics in all cases 
gives the value z — 3/2 corresponding to model F. This would then be different from the 3D 
XY model with spin dynamics in rcf. and the vortex loop prediction in r of. § . The other 
possibility is that all cases correspond to a vortex loop dynamics with a z slightly lower than 
3/2. Our present precision is not enough to distinguish between these alternatives. 

In conclusion we have from simulations determined the dynamic critical exponent z by using 
size scaling at T c for the 3D XY model with relaxational and RSJ dynamics. For relaxational 
dynamics with PBC we obtain z « 2 which is consistent with model A dynamics |t], fl2|| . 
However we conclude that the relaxational dynamics with FTBC, the lattice vortex loop 
model with relaxational dynamics and PBC|j), as well as RSJ dynamics with both PBC and 
FTBC all have the value z « 1.5. We suggest that the reason for this agreement is that all 
these models effectively corresponds to relaxational dynamics of the vortex loops. Model F 
corresponds to z = 3/2 which is consistent with our result although the slightly lower value 
z = 1.44 for vortex loops in ref. || is also consistent. 

This work was supported by the Swedish Natural Research Council through contract FU 
04040-332. 
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